Large-eddy simulation of a particle-laden turbulent channel flow

Large-eddy simulations of a vertical turbulent channel flow with 420,000 solid particles are performed in order to get insight into fundamental aspects of a riser flow The question is addressed whether collisions between particles are important for the ow statistics. The turbulent channel ow corresponds to a particle volume fraction of 0.013 and a mass load ratio of 18, values that are relatively high compared to recent literature on large-eddy simulation of two-phase ows. In order to simulate this ow, we present a formulation of the equations for compressible ow in a porous medium including particle forces. These equations are solved with LES using a Taylor approximation of the dynamic subgrid-model. The results show that due to particle-uid interactions the boundary layer becomes thinner, leading to a higher skin-friction coefcient. Important effects of the particle collisions are also observed, on the mean uid prole, but even more o on particle properties. The collisions cause a less uniform particle concentration\ud and considerably atten the mean solids velocity prole

Joint PDF modelling of turbulent flow and dispersion in an urban street canyon

The joint probability density function (PDF) of turbulent velocity and concentration of a passive scalar in an urban street canyon is computed using a newly developed particle-in-cell Monte Carlo method. Compared to moment closures, the PDF methodology provides the full one-point one-time PDF of the underlying fields containing all higher moments and correlations. The small-scale mixing of the scalar released from a concentrated source at the street level is modelled by the interaction by exchange with the conditional mean (IECM) model, with a micro-mixing time scale designed for geometrically complex settings. The boundary layer along no-slip walls (building sides and tops) is fully resolved using an elliptic relaxation technique, which captures the high anisotropy and inhomogeneity of the Reynolds stress tensor in these regions. A less computationally intensive technique based on wall functions to represent boundary layers and its effect on the solution are also explored. The calculated statistics are compared to experimental data and large-eddy simulation. The present work can be considered as the first example of computation of the full joint PDF of velocity and a transported passive scalar in an urban setting. The methodology proves successful in providing high level statistical information on the turbulence and pollutant concentration fields in complex urban scenarios.Comment: Accepted in Boundary-Layer Meteorology, Feb. 19, 200

Separation between coherent and turbulent fluctuations. What can we learn from the Empirical Mode Decomposition?

The performances of a new data processing technique, namely the Empirical Mode Decomposition, are evaluated on a fully developed turbulent velocity signal perturbed by a numerical forcing which mimics a long-period flapping. First, we introduce a "resemblance" criterion to discriminate between the polluted and the unpolluted modes extracted from the perturbed velocity signal by means of the Empirical Mode Decomposition algorithm. A rejection procedure, playing, somehow, the role of a high-pass filter, is then designed in order to infer the original velocity signal from the perturbed one. The quality of this recovering procedure is extensively evaluated in the case of a "mono-component" perturbation (sine wave) by varying both the amplitude and the frequency of the perturbation. An excellent agreement between the recovered and the reference velocity signals is found, even though some discrepancies are observed when the perturbation frequency overlaps the frequency range corresponding to the energy-containing eddies as emphasized by both the energy spectrum and the structure functions. Finally, our recovering procedure is successfully performed on a time-dependent perturbation (linear chirp) covering a broad range of frequencies.Comment: 23 pages, 13 figures, submitted to Experiments in Fluid

Fluid Particle Accelerations in Fully Developed Turbulence

The motion of fluid particles as they are pushed along erratic trajectories by fluctuating pressure gradients is fundamental to transport and mixing in turbulence. It is essential in cloud formation and atmospheric transport, processes in stirred chemical reactors and combustion systems, and in the industrial production of nanoparticles. The perspective of particle trajectories has been used successfully to describe mixing and transport in turbulence, but issues of fundamental importance remain unresolved. One such issue is the Heisenberg-Yaglom prediction of fluid particle accelerations, based on the 1941 scaling theory of Kolmogorov (K41). Here we report acceleration measurements using a detector adapted from high-energy physics to track particles in a laboratory water flow at Reynolds numbers up to 63,000. We find that universal K41 scaling of the acceleration variance is attained at high Reynolds numbers. Our data show strong intermittency---particles are observed with accelerations of up to 1,500 times the acceleration of gravity (40 times the root mean square value). Finally, we find that accelerations manifest the anisotropy of the large scale flow at all Reynolds numbers studied.Comment: 7 pages, 4 figure

Numerical ABL Wind Tunnel Simulations with Direct Modeling of Roughness Elements through Immersed Boundary Condition Method

Reproduction of atmospheric boundary layer wind tunnel experiments by numerical simulation is achieved in this work by directly modeling with immersed boundary method the geometrical elements placed in the wind tunnel's floor to induce the desired characteristics to the boundary layer.The wind tunnel has a cross section of 2.2 m x 2.25 m, with an inlet region 14 m long and a working region 2 m long. Boundary layer development is shaped up with a series of cubical elements, 3 cm in side, placed in a regular staggered arrangement with a 15 cm spacement. Vortex induction, Standen spires type elements, of 13,4 cm height, and a wall of 31.5 cm height are placed at the inlet. This arrangement is used to reproduce a representative urban site boundary layer (figure 1).The numerical model is implemented on the basis of the open source modelcaffa3d.MBRi [Usera et al 2008], which uses a finite volume method over block structured grids, coupled with various LES approaches for turbulence modeling and parallelization through domain decomposition with MPI [Mendina et al 2013]. Simulations were setup with approximately 2 million cells per block, with a 26 block arrangement. The computational grid is horizontally uniform with a resolution of 1.04 cm x 1.04 cm and nonuniform in vertical direction with the grid points concentrated near the floor . The grid spacing is geometrically stretched away from the floor with a minimum value of 1mm. The time step is 0.1 second and the computation is distributed in 26 cores on the Cluster-FING infraestructure [www.fing.edu.uy/cluster]. The Immersed boundary method approach followed the work of [Liao et al 2009]. Numerical simulation results are compared to wind tunnel measurements for the mean velocity profiles (figure 2), rms profiles and spectrums, providing good overall agreement. We conclude that the Immersed Boundary Condition method is a promising approach to numerically reproduce ABL Boundary Layer development methods used in physical modeling.Agencia Nacional de Investigación e Innovació

Macroscopic effects of the spectral structure in turbulent flows

Two aspects of turbulent flows have been the subject of extensive, split research efforts: macroscopic properties, such as the frictional drag experienced by a flow past a wall, and the turbulent spectrum. The turbulent spectrum may be said to represent the fabric of a turbulent state; in practice it is a power law of exponent \alpha (the "spectral exponent") that gives the revolving velocity of a turbulent fluctuation (or "eddy") of size s as a function of s. The link, if any, between macroscopic properties and the turbulent spectrum remains missing. Might it be found by contrasting the frictional drag in flows with differing types of spectra? Here we perform unprecedented measurements of the frictional drag in soap-film flows, where the spectral exponent \alpha = 3 and compare the results with the frictional drag in pipe flows, where the spectral exponent \alpha = 5/3. For moderate values of the Reynolds number Re (a measure of the strength of the turbulence), we find that in soap-film flows the frictional drag scales as Re^{-1/2}, whereas in pipe flows the frictional drag scales as Re^{-1/4} . Each of these scalings may be predicted from the attendant value of \alpha by using a new theory, in which the frictional drag is explicitly linked to the turbulent spectrum. Our work indicates that in turbulence, as in continuous phase transitions, macroscopic properties are governed by the spectral structure of the fluctuations.Comment: 6 pages, 3 figure

The effectiveness of neuromuscular warm-up strategies, that require no additional equipment, for preventing lower limb injuries during sports participation: a systematic review

PMCID: PMC3408383The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1741-7015/10/75. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited

The rise of fully turbulent flow

Over a century of research into the origin of turbulence in wallbounded shear flows has resulted in a puzzling picture in which turbulence appears in a variety of different states competing with laminar background flow. At slightly higher speeds the situation changes distinctly and the entire flow is turbulent. Neither the origin of the different states encountered during transition, nor their front dynamics, let alone the transformation to full turbulence could be explained to date. Combining experiments, theory and computer simulations here we uncover the bifurcation scenario organising the route to fully turbulent pipe flow and explain the front dynamics of the different states encountered in the process. Key to resolving this problem is the interpretation of the flow as a bistable system with nonlinear propagation (advection) of turbulent fronts. These findings bridge the gap between our understanding of the onset of turbulence and fully turbulent flows.Comment: 31 pages, 9 figure